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Technical  Summary 


An  experimentally  benchmarked  FORTRAN  computer  code  has  been 
developed  for  the  VAX  computer.  This  simulation  code  treats  the 
dynamic  processes  occurring  in  the  interior  of  electric  gun 
systems.  It  pays  special  attention  to  the  processes  occurring  in 
the  boundary  layer  between  the  hot  plasma  flow  and  the  solid 
"wall"  of  the  gun  barrel  bore. 

This  code  has  been  used  to  study  the  physical  processes 
involved  in  the  wall  erosion  problem  in  electric  gun  systems. 
This  includes  heat  transport  into  and  through  the  very  thin 
boundary  layer  and  the  effect  of  this  heat  transmission  on  the 
wall  material. 

An  analytical  model  of  the  flow  and  transport  in  an 
electrothermal  gun  has  been  developed.  This  model  accurately 
predicts  many  of  the  features  of  the  heat  transport  to  the  wall. 
This  process  is  the  primary  cause  of  erosion  in  the  gun  bore.  A 
report  describing  this  model  is  presented  as  Appendix  A. 

The  code  has  been  benchmarked  against  laboratory  experiments 
at  GT-Devices.  In  these  experiments,  plasma  flows  were  produced 
which  approximated  both  the  propellant  in  an  electrothermal  gun 
and  the  armature  in  a  rail  gun.  These  flows  were  directed  past 
samples  of  gun  bore  material,  and  the  ablation  of  the  material 
was  measured.  Since  the  properties  of  the  plasma  flow  were 
known,  the  heat  transport  through  the  skin  layer  to  the  sample 
could  be  calculated.  The  comparison  benchmarks  the  calculations. 

In  the  case  of  the  Electrothermal  gun,  the  study  of  these 
processes  has  led  to  a  change  in  the  design  of  the  gun  which  has 
virtually  eliminated  the  wall  erosion  problem  for  many  projectile 
masses  and  velocities  of  military  interest.  The  required  change 
is  the  addition  of  a  "mixing  chamber"  between  the  capillary  which 
is  the  source  of  hot  plasma  and  the  barrel,  which  requires  a 
cooler  environment. 

The  understanding  which  has  been  developed  for  the  wall 
erosion  problem  has  been  extended  to  the  Electromagnetic  (rail) 
gun  geometry.  The  heat  sources  internal  to  a  rail  gun  lead  to 
higher  heat  loads  on  the  rail  gun  wall  for  projectile  parameters 
of  military  interest.  This  heat  load  is  delivered  both  to  the 
metal  rails  and  to  the  insulators  between  them.  The  effect  on 

these  components  is  to  cause  ablation.  ,*-■ - . - — - - 
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Contract  Summary 


(a)  The  statement  of  work  for  this  program  is: 

This  principal  desired  outcome  of  this  overall  project 
is  an  experimentally  benchmarked,  VAX-compatible,  FORTRAN 
(or  PASCAL)  computer  simulation  code  which  will  treat  in  a 
physically  validated  fashion  the  dynamic  processes  occurring 
in  the  boundary  layer  between  a  dense,  high  temperature 
plasma  and  a  solid  "wall.''  The  code  will  employ  a  number  of 
asymptotically  coupled,  one-dimensional  algorithms  to  treat 
both  the  fast  timescale  processes  perpendicular  to  the  wall 
as  well  as  the  much  slower  timescale  plasma  flow  along  the 
wall.  This  code  will  be  used  to  enhance  the  physical 
understanding  of  the  wall  erosion  problem  in  electrothermal 
guns  and  will  be  available  for  possible  application  to  the 
rail  erosion  problem  in  electromagnetic  launchers  by  other 
research  efforts.  The  desired  research  milestones  for  each 
of  the  three  years  include  the  following: 

Determine  the  appropriate  transport  coefficients  to  be 
used  in  the  electrothermal  jet  and  in  the  wall  boundary 
layer  including  the  Rosseland  opacity  for  radiation 
transport,  the  viscosity,  and  the  electrical  and  thermal 
conductivities  for  the  partially  ionized  plasma.  Also, 
integrate  these  transport  coefficients  into  a  VAX- 
compatible,  FORTRAN  (or  PASCAL)  code  which  contains  the  flow 
algorithms  for  the  boundary  in  the  gun  barrel.  The 
resultant  code  will  be  used  to  simulate  systems  with 
experimentally  relevant  parameters  and  the  results  will  be 
compared  to  empirical  data  available  from  electrothermal 
guns  presently  in  operation.  An  annual  report  will  be 
delivered  to  AFOSR  which  will  discuss  the  algorithms 
contained  in  the  code,  will  summarize  the  findings  of  the 
comparison  with  experimental  data  and  will  include  a 
complete  copy  of  the  FORTRAN  (or  PASCAL)  source  file  of  the 
code(s)  on  an  IBM-PC  compatible  diskette. 

Option  1: 

Based  on  the  experimental  comparisons,  evaluate  the 
importance  of  turbulence  in  the  erosion  process.  If 
necessary,  incorporate  a  turbulent  boundary  layer  algorithm 
in  the  code.  Conduct  a  series  of  simulations  for  a  range  of 
electrothermal  gun  parameters,  electrical  power  levels,  and 
materials.  Determine  the  dependencies  of  the  gun 
performance  parameters  to  these  changes.  Submit  an  annual 
report  to  AFOSR  to  present  the  major  results  of  the  above 
work.  The  report  should  include  an  IBM-PC  compatible 
diskette (s)  containing  a  complete  copy  of  the  FORTRAN  (or 
PASCAL)  source  file  of  the  latest,  updated  version  of  the 
simulation  code(s). 
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Option  2: 

Incorporate  into  the  code(s)  output  routines  which  will 
generate  graphical  displays  of  the  solutions  so  that  the 
dynamics  can  be  visualized.  Use  the  computer  simulation 
results  to  derive  some  general  scaling  relationships  for  the 
dynamic  behavior  of  the  plasma  jets.  Submit  a  final  report 
for  AFOSR  which  will  include  a  compilation  of  simulation 
results  and  the  derived  scaling  relationships  that  would  be 
useful  for  applications  of  these  electrothermal  jets  to 
hypervelocity  launchers  and  to  other  applications.  This 
final  report  will  include  a  complete  copy  of  the  FORTRAN'  (or 
PASCAL)  source  file  of  the  latest,  updated  version  of  the 
simulation  code  on  an  IBM-PC  compatible  diskette  as  well  as 
complete  documentation  for  the  code's  algorithms  and 
instructions  for  its  use. 

(b)  Status  of  the  research  effort: 

All  the  milestones  of  the  initial  task  and  option  1 
have  been  met.  Option  2  has  not  been  exercised.  Activity 
under  the  contract  has  been  completed.  The  experimentally 
benchmarked,  VAX-compatible  FORTRAN  computer  simulation  code 
has  been  developed.  This  code  was  used  to  enhance  the 
physical  understanding  of  wall  erosion.  The  transport 
coefficients  for  Rosseland  opacity,  viscosity,  electrical 
and  thermal  conductivity  and  turbulence  have  been 
incorporated.  Simulations  have  been  performed  which  are 
relevant  to  both  electrothermal  guns  and  electromagnetic 
(rail)  guns.  This  final  report  has  been  submitted.  A 
complete  copy  of  the  FORTRAN  source  file  is  included  with 
the  copies  of  the  report,  on  a  DOS  3  formatted  5-1/4"  360k 
diskette . 

(c)  No  journal  articles  have  been  prepared  or  published;  however, 

a  technical  report, 

"Gun  Interior  Ballistics  Model  for  Constant  Acceleration 
Launch  of  Projectiles,"  by  D.  A.  Tidman,  N.  K.  Winsor  and 
S.  A.  Goldstein,  GT-Devices  Technical  Note  GTD  86-2,  dated 
February  1986, 

has  been  prepared  and  distributed.  It  is  included  here  as 

Appendix  A. 

(d)  The  professional  personnel  associated  with  the  program  are: 

Dr.  Niels  K.  Winsor 
Dr.  Derek  A.  Tidman 
Dr.  Shyke  A.  Goldstein 
Dr.  Howard  W.  Bloomberg 
No  advanced  degrees  have  been  awarded. 
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(e) 


Interactions : 

1.  Paper  presented  by  N.  K.  Winsor,  titled  "The  Gun  Bore 
Bearing  in  EM  and  ET  Launchers,"  by  N.K.  Winsor,  D.A. 
Tidman.S.A.  Goldstein  and  H.W.  Bloomberg,  presented  at 
Innovative  Science  and  Technology  Electromagnetic 
Launch  Science  Interchange  Workshop,  Eglin  AFB,  23  June 
1986. 

2.  Paper  presented  by  D.  A.  Tidman,  titled  "A  Rail  Gun 
Plasma  Armature  Model , "  at  the  Electromagnetic  Gun 
Armature  Workshop,  Eglin  Air  Force  Base,  24-26  June 
1986. 

3.  Half-day  seminar-workshop  presented  by  N.  K.  Winsor  and 
Y.  Chia  Thio,  to  Eglin  AFB  and  Auburn  University 
personnel  on  "Electromagnetic  and  Electrothermal 
Research  at  GT-Devices,"  Eglin  Air  Force  Base,  21-22 
August  1986. 


(f)  New  discoveries,  inventions,  or  patent  disclosures  and 
specific  applications  stemming  from  the  research  effort: 

In  conjunction  with  experiments,  these  simulations  aided  in 
interpreting  the  erosion  phenomena  ih  the  GT-Devices 
9  20  mm  Electrothermal  gun.  The  synergism  of  these  studies 

provided  the  insight  necessary  to  redesign  the  breech 
chamber  of  the  gun  and  substantially  eliminate  the  erosion 
problem  for  most  applications  of  military  interest. 


(g) 


Any  other  statements: 
See  balance  of  report! 
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Introduction 


The  performance  of  an  airborne  gun  significantly 
affects  the  kind  of  mission  the  aircraft  can  accomplish.  For 
example,  the  C-130  gunship  is  designed  for  use  in  an  orbit  around 
a  selected  target.  A  higher  velocity  gun  can  mean  the  difference 
between  an  orbit  within  range  of  ground  fire  and  a  position  well 
outside  the  range  of  the  ground  defense.  (See  Figure  1.)  The 
higher  velocity  can  also  simplify  fire  control  when  engaging  a 
moving  target. 


For  conventional  powder  guns,  a  practical  velocity  limit  for 
long-life  service  is  around  1  km/sec  (3300  ft/sec).  For  example, 
this  is  close  to  the  performance  of  the  Phalanx  close-in 
shipboard  defense  system,  which  has  been  designed  for  high 
velocity  and  rapid  fire. 


Chemical  propellants  are  capable  of  higher  velocities,  but 
when  such  performance  is  achieved,  it  is  at  considerable  cost, 
usually  in  a  much  larger  chemical  charge,  greatly  increased 
barrel  wear  and  complex  projectile  designs,  often  including  a 
discarding  sabot. 


The  Strategic  Defense  Initiative  has  motivated  the  search 
for  methods  to  achieve  much  higher  velocity  with  gun  systems. 
One  result  has  been  increased  research  on  electric  gun  systems, 
principally  the  electromagnetic  (rail)  gun. 


More  recently,  the  electrothermal  gun  has  been  receiving 
attention  as  a  gun  with  many  of  the  advantages  of  the  rail  gun, 
while  avoiding  most  of  the  disadvantages  of  other  kinds  of 
electric  guns. 
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Electric  guns  share  advantages  such  as  higher  velocities, 
insensitive  munitions,  velocity  control  from  round  to  round  and 
more  rounds  storable  in  a  given  magazine  size. 

The  disadvantages  come  primarily  in  the  requirement  for  an 
electrical  power  conditioning  system  and  in  increased  complexity 
of  the  complete  system.  If  an  electrical  power  supply  is  already 
available,  an  electric  gun  may  be  a  convenient  add-on  to  other 
systems . 

There  are  significant  differences  between  the  two  main 
electric  gun  technologies.  Electrothermal  gun  can  use  standard 
gun  components  and  munitions,  while  rail  guns  require  a  split 
barrel  and  nonconducting  (or  insulated)  munitions.  On  the  other 
side,  rail  guns  are  theoretically  able  to  achieve  higher 
velocities,  because  the  driving  force  follows  the  projectile  down 
the  barrel. 

Both  kinds  of  electric  guns  have  potential  problems  with 
barrel  erosion.  Control  of  this  erosion  will  improve  the 
accuracy  of  the  gun  and  the  lifetime  of  the  components. 

The  objective  of  this  effort  is  to  produce  analytic  and 
numerical  models  of  the  processes  which  are  related  to  this 
erosion  and  to  the  performance  of  the  guns.  The  gas  flow  along 
the  barrel  wall  and  the  heat  transport  to  that  wall  play  a 
crucial  role  in  these  processes. 

The  following  chapters  describe  in  turn  the  physics  of  the 
flow  and  transport  processes,  their  numerical  treatment  and  the 
results  of  computations.  These  are  followed  by  a  chapter  drawing 
conclusions  from  this  year  of  work. 


Physics 


Conventional  service  guns  have  relatively  low  velocity 
capability,  but  specialized  guns  have  no  such  limit.  So-called 
"light-gas  guns"  have  achieved  velocities  over  5  km/sec  with 
kilogram  projectiles,  and  approach  10  km/sec  with  few-gram 
projectiles.  They  are  not  weaponizable  because  some  components 
of  the  breech  are  destroyed  in  each  shot. 

These  light-gas  guns  provide  considerable  information  about 
the  limits  of  component  performance  in  high-velocity  guns.  The 
gas  flow  in  the  barrel,  and  the  damage  it  does  to  the  barrel  are 
not  well  understood. 

The  basic  flow  process  is  straightforward.  A  source  of  mass 
and  energy  at  the  breech  end  of  the  gun  produces  high  pressure. 
This  pressure  accelerates  the  propellant  material  down  the 
barrel,  accelerating  a  projectile. 

Two  things  limit  the  ultimate  speed  this  projectile  can 
reach.  One  is  the  sound  speed  of  the  gas.  The  second  is  the 
interaction  between  the  gas  flow  and  the  barrel  wall. 


Gradual  or  Abrupt 
Reacted  Propellant  '  Transition  Section 


Figure  2.  Schematic  illustration  of  conventional  gun  structure. 
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The  effect  of  the  gas  sound  speed  is  well  known  and  readily 
calculable.  Figure  2  shows  the  basic  structure  of  a  conventional 
gun.  Reacted  propellant  of  pressure  P  and  temperature  T  in  a 
chamber  reservoir  expand  through  a  breech  and  accelerate  a 
projectile  of  mass  M  down  a  barrel.  This  process  bears  a  close 
relation  to  the  operation  of  an  Electrothermal  gun,  and  a  more 
distant  relation  to  the  operation  of  an  Electromagnetic  gun. 

Seigel1  has  calculated  the  performance  of  a  "Pre-burned 
Propellant  Ideal  Gas"  (PPIG)  gun,  and  numerically  determined  the 
acceleration  that  can  be  expected  with  various  gun 
configurations.  His  results  are  in  excellent  agreement  with 
experiments  on  chemically  driven  high-velocity  guns  such  as 
light-gas  guns. 

His  analysis  begins  with  a  simplified  version  of  the 
equations  which  are  presented  in  the  following  section.  This  is 
a  natural  starting  place  for  the  analysis  of  any  high-velocity 
gun.  These  equations  will  be  used  to  analyze  the  dynamics  of 
electric  guns. 


General  equations 


The  propellant  in  a  conventional  gun  is  obtained  as  the 
reaction  products  of  a  chemical  burn  (the  PPIG  pre-burned 
propellant) .  In  an  electrothermal  gun,  the  propellant  is 
evaporated  from  fill  material  in  a  capillary  tube  behind  the  gun 
breech  (Figure  3) .  This  means  the  pressure  and  temperature  of 
the  propellant  are  determined  by  the  electrical  energy  input, 
rather  than  simple  adiabatic  expansion. 
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Figure  3.  Electrothermal  gun  capillary  structure. 

In  an  electromagnetic  gun,  the  propellant  is  composed  of  armature 
and  barrel  material  ablated  by  the  armature  arc.  The  electrical 
energy  is  deposited  near  the  projectile  rather  than  at  the 
breech,  but  the  force  is  still  transmitted  to  the  projectile  by 
propellant  pressure. 


In  general,  gun  material  is  converted  to  propelling  fluid 
and  compressed  by  thermal  or  magnetic  forces.  The  propelling 
fluid  satisfies  the  continuity  equation 


P  ¥ 


(3.1) 


where  6p/5t  is  the  rate  of  ablation  or  evaporation  of  propelling 
fluid.  In  the  barrel,  the  wall  ablates,  creating  the  erosion 
problem  that  limits  barrel  lifetime.  In  addition  to  fluid 
production,  pressure  gradients  and  viscous  stress  affect  the 
fluid  flow  velocity 

p^.JxB-£P-V-I  ,3.2) 

and,  in  conjunction  with  other  transport  processes,  they 
determine  the  evolution  of  the  fluid  total  energy 
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{if  ♦  E  2  •  v  .  !£  -  2  •  v)  +  W  -  Q 

where 

E  =  i  PV2  +  U  +  P 

and  U  is  the  fluid  internal  energy.  The  heat  input  is  by  Ohmic 
heating  where  current  flows  in  the  fluid 


W  =  n  J2 


and  the  heat  transfer  and  loss  terms  considered  here  are  thermal 
conduction  and  radiation 

Q  =  -  V  •  I^VT  -  V  •  KR VT 

and  turbulence.  The  turbulence  has  been  examined  in  detail 
analytically.2  These  studies  have  shown  that  the  effects  of 
small  eddies  can  be  approximated  by  a  term  with  the  same 
functional  dependence  as  the  bulk  viscosity.  The  result  of  such 
an  inclusion  is  the  production  of  larger-scale  eddies  which  are 
very  similar  to  those  which  are  produced  by  conventional 
turbulent  mixing  in  the  boundary  layer. 

These  equations  describe  the  dynamics  of  the  propelling 
fluid  flow,  but  they  must  be  supplemented  by  information  about 
the  properties  of  the  fluid  itself.  This  information  is  provided 
by  the  equation  of  state.  We  will  use  it  in  the  form 

T  =  T(p,U)  (3.4) 

P  =  P(p,U) 

After  specifying  these  quantities,  the  mathematical  problem  is 
complete  when  the  transport  coefficients  are  specified. 
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The  fluid  heating,  W,  in  both  the  ET  and  EM  case  is  by 
electrical  Ohmic  heating.  The  plasma  has  a  resistivity3 

n  =  nQ  Z£nA(l  +  2  vg0  /ve-  )/nJ3/2  (3.5) 

The  location  of  the  heating  is  different  in  the  two  cases  and 
will  be  discussed  in  the  following  sections. 


The  viscosity  is  a  symmetric  tensor  proportional  to  the 
velocity  shear  with  geometric  adjustments  for  divergent  flow4 


w  -  tc  r»Va 
[axe 


ave 

3Xa 


2 

3  6ae 


V  •  vj 


(3.6) 


with  a  coefficient  that  in  weak  magnetic 
is  a  strong  function  of  temperature 
density3 


5/2 

|f  =  0.42  (K-T) - 

^  e4  £n  A 


fields  or  high  density 
and  a  weak  function  of 


(3.7) 


Its  effect  on  the  energetics  of  the  flow  is  to  convert  kinetic 
energy  to  thermal  energy  in  regions  of  sheared  flow. 


In  regions  with  substantial  Ohmic  heating,  the  central  bore 
temperature  goes  well  above  the  melting  point  of  the  wall,  and 
the  boundary  layer  is  subject  to  heat  transport  by  both 
conduction  and  radiation.  The  radiation  transport  is  a  complex 
non-local  process,  but  in  the  black-body  case  it  can  be  modeled 
by  a  diffusion  coefficient 


qR  =  -  KRVT 


(3.8) 
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where  Kg  is  the  Rosseland  mean  opacity.  It  is  tabulated  in  the 
Los  Alamos  National  Laboratories  SESAME  Tables.8 


Electrothermal  case 

An  electrothermal  gun  is  a  natural  extension  of  classical 
gun  technology.  It  uses  conventional  gun  barrels  and  projectiles. 
The  new  advance  comes  from  replacing  the  chemical  propellant  with 
a  combination  of  a  propelling  medium  of  choice  and  electrical 
energy . 

Figure  4  is  an  illustration  of  an  electrothermal  (ET)  gun. 
Electric  energy  is  delivered  to  the  gun  by  a  power  conditioning 
system  that  shapes  the  electrical  power  pulse  as  needed  for  gun 
operation.  A  cartridge  confines  an  electrical  arc  which  converts 
the  power  into  heat,  pressurizing  the  propelling  fluid. 

After  the  electrical  energy  has  been  delivered  to  it,  the 
gun  behaves  like  a  conventional  powder  gun,  with  the  propellant 
chasing  the  projectile  down  a  conventional  metal  barrel.  The  ET 
gun  greatly  outperforms  a  chemical  gun  for  two  reasons: 

1.  The  propelling  fluid  may  be  chosen  to  have  better  fluid 
properties  than  chemical  propellants,  since  it  is  not 
restricted  to  exothermal  materials.  In  particular,  it  can 
be  a  low  molecular  weight  material,  to  produce  a  high  sound 
speed  and  therefore  a  propelling  fluid  which  can  transmit 
pressure  down  the  barrel  at  high  speeds. 
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Figure  4.  Illustration  of  Electrothermal  gun  structure. 


2.  The  electrical  power  pulse  can  be  shaped  to  maintain  high 
pressure  throughout  the  projectile  acceleration  down  the 
barrel.  This  is  in  contrast  to  chemical  burn,  which  usually 
produces  a  rapid  rise  to  peak  pressure  before  the  projectile 
has  had  much  opportunity  to  move,  followed  by  a  decaying 
pressure  as  the  projectile  accelerates.  The  ET  power  input 
can  be  programmed  to  produce  a  constant  pressure  near  the 
maximum  pressure  the  gun  can  withstand. 

Figure  5  compares  the  theoretical  and  measured  performance  of  the 
ET  gun  with  conventional  guns.  The  ET  gun  can  reach  speeds 
unavailable  to  conventional  guns ,  and  do  reach  comparable 
velocities  with  much  higher  efficiency. 

The  ET  gun  also  has  problems  which  must  be  considered.  The 
principal  one  is  that  the  high  temperatures  in  the  cartridge,  and 
the  high  flow  velocities  can  easily  erode  and  damage  the  gun 
bore.  Powder  guns  also  have  a  bore  erosion  problem,  but  they 
also  have  many  years  of  design  and  development  devoted  to  its 
management.  The  present  effort  has  aided  similar  developments 
for  the  electrothermal  gun. 

The  energy  is  delivered  to  the  ET  gun  by  Ohmic  heating  in 
the  cartridge.  In  this  region,  the  temperatures  reach  typically 
40,000° K,  and  the  surface  of  solids  rapidly  evaporates.  The 
densities  are  high  enough  that  the  radiation  takes  on  black-body 
form  and  the  central  temperature  becomes 

T(eV)  =  1.8 

This  relation  is  obtained  from  steady-state  equations  for  the 
heat  deposition  and  transport  in  the  cartridge. 
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Figure  6  schematically  indicates  the  shape  of  this  layer, 
and  the  competing  transport  processes  which  determine  the  shape 
and  width  of  the  layer.  Because  of  the  very  steep  gradients  in 
this  sheath,  the  transport  to  the  wall  takes  place  on  a  much 
shorter  time  scale  than  the  axial  transport.  This  allows  the 
radial  transport  to  be  treated  asymptotically  in  the  axial 
transport  calculations. 


Figure  6.  Thermal  transport  skin  layer  in  a  gun  bore. 


Electromagnetic  case 


An  electromagnetic  (EM)  gun  is  depicted  schematically  in 
Figure  7.  Like  the  ET  gun,  the  electric  energy  is  delivered  at 
the  breech  of  the  gun  (or,  in  more  complicated  arrangements,  at  a 
number  of  feed  points  along  the  bore).  Unlike  the  ET  gun,  the 
EM  gun  deposits  electrical  energy  close  to  the  projectile,  with 
the  deposition  region  traveling  along  the  bore  with  it. 

The  bore  is  split,  consisting  of  two  rails  and  two 
insulators,  usually  arranged  in  a  circular  or  square  cross- 
section.  It  is  normally  under  external  compression,  at  least  to 
the  peak  pressure  anticipated  during  gun  operation. 

The  energy  deposition  region  within  the  bore  is  called  an 
"armature"  in  analogy  to  the  armature  on  a  motor,  where  current 
is  transferred  from  a  stationary  conductor  to  a  moving  conductor. 
In  the  1970 's,  solid  armatures  were  tried,  but  it  was  found  that 
a  thin  plasma  formed  between  the  solid  and  the  rails,  and  damaged 
the  rails.  The  damage  could  be  reduced  by  increasing  the  size 
(and  mass)  of  the  armature,  but  that  reduced  the  useful 
projectile  payload.  As  a  result,  most  high-velocity  railguns  now 
use  plasma  armatures,  both  to  reduce  the  armature  mass  and  to 
increase  the  area  of  the  armature,  thus  hopefully  spreading  out 
the  heat  load  on  the  bore. 

In  the  EM  gun,  the  current  flows  from  one  rail  to  the  other. 
In  this  geometry,  the  resistance  is 

R  =  nQ  Z  tnA  (1  +  2  veo/vei)  G/xAT3/2  (3.10) 

where  xa  is  the  armature  length  and  the  geometry  is  characterized 
by  the  factor  G  =1  for  a  square  bore,  and  G  =  n/4  for  circular 
bore.  For  a  constant-mass  armature  with  sound  speed  cs  ,  the 
armature  length  is 


CURRENT 


Figure  7.  Structure  of  an  Electromagnetic  (rail)  gun. 
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xa  =  n  V 

where  A  is  the  cross-sectional  area  of  the  bore,  and  )'  and  Ma  are 
the  armature  adiabatic  constant  and  mass,  respectively. 

The  propelling  fluid  in  the  EM  gun  may  start  out  as  some 
fuse  or  other  selected  material,  but  it  usually  quickly  becomes 
dominated  by  rail  and  insulator  material.  The  bore  wall  is 
subjected  to  the  very  high  temperatures  of  the  armature  plasma. 
As  the  armature  moves  down  the  bore,  it  evaporates  wall  material, 
which  may  be  swept  up  in  the  armature  arc  and  carried  along  with 
the  projectile.  This  mix  of  materials  is  the  armature  fluid. 

The  primary  focus  of  this  research  effort  is  the  flow  and 
operation  of  the  electrothermal  gun.  However,  the  models  which 
have  been  developed  retain  sufficient  generality  to  treat  the 
electromagnetic  gun  with  appropriate  changes  in  the  resistivity 
and  geometry. 
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Numerics 


The  equations  presented  in  the  previous  section  are  both 
complex  and  nonlinear.  Some  progress  can  be  made  in  their 
solution  analytically,  but  many  approximations  are  required. 

Numerical  calculation  of  the  solutions  has  the  advantage 
that  accurate  results  can  be  obtained  without  these 
approximations.  It  has  the  disadvantage  that,  for  nonlinear 
equations,  a  numerical  solution  applies  only  for  the  specific 
case  calculated,  and  not  for  a  class  of  solutions,  as  can  be 
obtained  with  analytic  scaling  laws. 

On  the  other  hand,  numerical  results  can  often  lead  to 
empirical  relations  which  are  very  useful,  and  can  verify  the 
accuracy  of  analytic  results  which  are  obtained  via 
approximations  whose  validity  is  otherwise  hard  to  assess. 

This  section  takes  the  equations  from  the  previous  section, 
converts  them  to  forms  which  are  more  conveniently  handled 
numerically,  and  sketches  how  they  are  solved  numerically.  This 
provides  the  framework  within  which  the  transport  research  and 
electric  gun  development  is  carried  out. 


Conservation  laws 

Numerical  simulation  of  fluid  flow  in  the  gun  bore  has  the 
advantage  that  it  can  treat  non-linearities  and  functional 
dependencies  without  approximation.  It  has  the  disadvantage  that 
it  makes  other  approximations  to  treat  a  problem  with  practically 
infinite  degrees  of  freedom.  One  way  or  another,  it  reduces  the 
problem  to  a  finite  number  of  degrees  of  freedom. 
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The  method  used  here  is  finite-difference  formulation  of  the 
numerical  problem.  The  continuous  variables  of  the  problem  are 
replaced  by  a  finite  number  of  discrete  values,  arrayed  along 
coordinate  axes.  It  is  then  desirable  to  express  the  equations 
in  a  way  which  preserves  the  physics  of  the  problem  as  well  as 
possible,  given  the  necessity  for  this  discretization. 

One  effective  method  is  to  cast  the  equations  to  be  solved 
in  conservation-law  form.  This  is  easy  for  the  mass,  momentum 
and  energy  equations.  Then  a  numerical  method  can  be  used  which 
preserves  this  conservation  property  to  the  maximum  accuracy 
available  on  the  computer  of  choice. 

This  is  how  we  proceed  here. 

The  mass  transport  equation  (3.1)  can  be  converted  to 
conservation-law  form  by  rewriting  its  total  time  derivative  as  a 
partial  time  derivative  plus  a  convective  derivative  and 
collecting  terms 

I|  +  H-pv=.|£  <4.i> 

Multiplying  this  equation  by  v  and  subtracting  it  from  Eq.  (3.2) 
yields  the  conservation-law  form  of  the  momentum  equation 

"  (pv)  +  V  •  (pyv)  =  J  x  B  -  VP  -  V  •  n  (4.2) 

A  similar  manipulation  converts  the  energy  equation  to 

fl  *  -  '  Eu  -  -  v-  («-v)  +  W  -  Q  (4'3) 

The  remainder  of  the  equations  in  Section  3  are  scalar,  and  do 
not  require  further  manipulation. 
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Finite-difference  formulation 
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The  simplest  view  of  the  finite  difference  formulation  of 
these  differential  equations  can  be  illustrated  by  expressing  the 
derivatives  in  Eq.  (4.1)  in  terms  of  values  of  the  variable  on  a 
discrete  one-dimensional  mesh.  That  is. 


V-  PV  =  (pH1 


vi+l 


+  0  (X1+l  -  Xi)  2 


(4.4) 


where  the  subscript  i  refers  to  the  value  of  the  variable  at  the 
location  xi  .  By  the  mean  value  theorem,  this  is  the  exact  value 
of  the  derivative  at  some  point  between  xi  and  xi  +  i  .  The  most 
probable  location  is  the  midpoint 
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i+i 


(4.5) 


The  next  step  in  a  finite-difference  formulation  is  discretizing 
the  time.  In  terms  of  a  time-step  6t,  the  mass  conservation 
equation  becomes 


(p'<‘+{t)-  p<l>)  /  6t  ♦  (pw  Vm 


6p 

6t 


(4.6) 


Again,  the  "centering"  of  the  time  values  is  important.  With 
centered  differences  in  both  time  and  space,  this  equation  is 


/  (t+fit) 
\i+* 


fit  + 


(  (t+i<5t)  (t+|5t)  (t+i6t)  \t+i6t) 

P  V  -  p  V 

i  i+1  i+1  1  i 


X  )  = 

A4  /  At 


fit  (4.7) 


This  result  is  the  simplest  second-order  accurate  finite- 
difference  formulation  of  Eq.  (3.1).  It  is  called  a  "leap-frog" 
algorithm  because  the  density  leaps  from  time  t-5t/2  to  t+5t/2 
over  the  values  of  density  and  velocity  at  time  t. 

One  major  problem  with  using  a  formulation  such  as  this  is 
that  the  values  of  density  and  velocity  must  be  determined  at 
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different  times  and  locations,  and  unless  great  care  is  taken  at 
keeping  them  synchronized,  these  different  grids  can  drift  apart, 
yielding  nonphysical  results. 

For  this  reason  it  is  prudent  to  use  a  numerical  algorithm 
which  has  been  thoroughly  tested  and  designed  for  the  kind  of 
solutions  one  expects  from  the  equations.  The  Flux-Corrected 
Transport  algorithm  is  chosen  for  this  case. 


Flux-corrected  transport 

The  propelling  fluid  in  the  gun  bore  has  a  number  of 
properties  that  can  make  numerical  treatment  difficult.  It  has 
steep  gradients  near  the  boundaries.  The  electrical  discharge 
can  produce  rapid  time  rates  of  change,  causing  dynamic  steep 
gradients,  perhaps  leading  to  shock  waves. 

Flux-Corrected  Transport  is  an  algorithm  which  was  developed 
with  such  applications  in  mind.  Figure  8  compares  it  with 
leapfrog  and  two  other  popular  algorithms  in  the  propagation  of  a 
very  steep  (square)  density  wave. 

FCT  begins  with  a  second-order  accurate  time-advance 
algorithm,  and  extends  it  by  examining  the  flux  at  each  grid 
point  to  insure  that  certain  physically-derived  inequalities  are 
satisfied.  (These  inequalities  are  violated  by  many  popular 
numerical  algorithms  when  steep  gradients  or  shock  behavior  is 
obtained  in  the  solution.  The  result  of  the  violation  is  non¬ 
physical  oscillatory  behavior  of  the  numerical  solutions.) 

The  FCT  algorithm  was  published  in  the  form  of  a  FORTRAN 
program.  It  is  used  here  with  minor  variations  on  the  published 
form . 
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The  FCT  algorithm  assumes  that  the  equations  have  been  cast 
in  conservation-law  form,  and  calculates  the  time-advance  of  each 
conserved  quantity,  given  the  right-hand  side,  or  "source"  terms. 
That  is,  FCT  is  formulated  principally  as  a  Fortran  subroutine 
which  accepts  as  input  vectors  of  conserved  quantities  and  source 
terms,  and  produces  as  output  the  vector  at  a  future  time. 

One  very  important  feature  of  FCT  is  that  it  makes  use  of 
all  quantities,  both  variables  and  sources,  at  the  same  time  and 
space  positions,  avoiding  the  problem  described  above  in  terms  of 
the  leapfrog  algorithm. 

In  order  to  obtain  optimum  performance  from  FCT,  and  to 
preserve  second-order  accuracy  in  time,  it  is  used  in  a 
predictor-corrector  formulation.  That  is,  FCT  calls  are  first 
made  to  advance  each  of  Eqs .  (4. 1-4. 3)  a  time  &6t.  Then  the 
sources  are  evaluated  and  a  second  series  of  calls  is  made  to 
advance  the  physical  variables  to  the  full  time  step,  6t. 
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Computation  and  experiment 


For  electric  guns,  both  modeling  and  experiment  have 
advanced  substantially  during  this  year.  Modi f ications  of 
Electrothermal  guns  have  allowed  operation  with  virtually 
ablation-free  barrels.  In  Electromagnetic  guns,  the  thermal 
behavior  of  the  conducting  rails  and  the  insulators  has  been 
characterized,  leading  to  a  more  complete  understanding  of  how  to 
operate  within  the  limitations  of  the  materials. 

The  computer  model  has  been  written  in  a  general  form  to 
allow  calculations  for  both  kinds  of  guns.  It  is  organized  with 
distinct  physical  processes  in  separate  subroutines.  Each 
subroutine  has  a  specific  function,  divided  by  physical  or 
computational  processes,  such  as  electrical  heating,  fluid 
properties,  acceleration,  transport,  diagnostics  and  output. 

Each  subroutine  has  multiple  entries.  The  different  entries 
in  each  subroutine  perform  parallel  functions,  such  as 
initialization,  data  entry,  parameter  evaluation,  output  and 
termination.  Compared  with  the  use  of  subroutines  without 
entries,  the  use  of  multiple  entries  reduces  the  use  of  COMMON 
statements  for  the  variables. 

Nevertheless,  most  of  the  key  variables  are  stored  in  COMMON 
blocks.  These  blocks  are  named,  and  entered  into  the  subroutines 
that  use  them  by  INCLUDE  statements.  This  prevents  misspellings 
or  similar  problems  which  could  put  the  COMMON  variables  at 
different  locations  in  different  subroutines. 
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Program  structure 


The  program  is  written  in  FORTRAN  for  the  Digital  Equipment 
Corporation  VAX  computer  series.  Because  the  VAX  is  a 
minicomputer,  some  compromises  must  be  made  in  the  program 
design . 

For  convenient  access  to  the  source,  it  is  desirable  to  have 
the  program  in  a  single  file.  This  allows  such  convenient 
operations  as  global  search,  global  replace,  block  moves  and 
other  editing  practices  which  are  more  conveniently  done  on  a 
single  file. 

On  the  other  hand,  a  model  of  this  size  takes  a  long  time  to 
develop,  and  a  long  time  to  completely  recompile  after  a  change 
has  been  made.  For  this  reason,  separate  modules  might  be 
preferred,  so  that  when  only  one  or  a  few  have  been  changed,  only 
those  few  must  be  recompiled  and  linked. 

To  make  both  of  these  approaches  possible,  the  program  was 
structured  so  that  both  options  were  available.  A  command  file 
was  built  containing  the  source  routines,  separated  by  CREATE 
commands.  The  created  files  have  alphabetically  ordered  names  so 
that  the  contained  subroutines  can  be  found  in  the  same  order  in 
the  VAX  directory. 

When  complete  recompilation  is  desired,  the  command  file  is 
executed,  splitting  the  complete  source  into  a  number  of  files, 
compiling  all  of  them  and  storing  the  object  files  in  a  library. 
The  smaller  source  files  are  then  available  for  change  and  rapid 
recompilation . 

In  this  way,  the  complete  executable  file  can  be  built  from 
individually  recompiled  files  and  the  object  modules  which  are 
stored  in  the  library.  This  speeds  development  when  small  parts 
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of  the  program  are  changed  and  global  operations  are  not 
required . 


Electrothermal  results 


Heat  flow  to  the  wall  of  a  gun  barrel  has  been  studied 
empirically  for  many  years.  More  recently,  progress  has  been 
made  in  analytically  determining  the  importance  of  various 
factors  in  the  process.6 

For  the  Electrothermal  gun,  one  main  objective  of  heat 
transport  analysis  is  to  identify  the  conditions  under  which  the 
gun  can  be  operated  without  erosion  of  the  gun  bore.  The  thermal 
skin  in  which  all  the  heat  transport  processes  are  competing  is 
the  key  to  understanding  the  physics  of  this  process. 

Figure  6  illustrates  the  interplay  of  the  major  processes  in 
the  transport.  Radiation  carries  energy  toward  the  wall.  As 
transport  to  the  wall  begins  to  cool  the  edge  of  the  propelling 
gas  flow,  the  black-body  transport  gives  way  to  radiation 
diffusion,  and  finally  direct  radiation  to  the  wall.  Radiation 
is  competing  in  the  thermal  sublayer  with  direct  thermal 
conduction  and  turbulent  convection,  both  of  which  become 
progressively  more  rapid  as  the  temperature  gradient  becomes 
steeper  at  the  bore  wall. 


The  response  of  the  wall  to  the  heat  load  from  the  flowing 
propellant  is  described  by  a  heat  flow  equation  for  the  metal  of 
the  barrel.  It  is 


r  9'T 
P  t  — 

w  v  at 


K-r  2  “  £l, 
3x 


w 


(5.1) 
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where  T(x,t)  is  the  wall  surface  temperature ,  Pw  is  the  wall  mass 
density,  Cv  is  its  heat  capacity,  qw  is  the  heat  load  on  the  wall 
and  Kv  is  its  thermal  conductivity.  The  solution  of  this 
equation  requires  advanced  methods,  such  as  Green's  functions7  or 
asymptotic  expansion  of  the  Riemann  function.8  The  key  result  is 
that  the  homogeneous  solution  ( for  qw  =  0)  of  Eq.  (5.1)  is 


TH  = 


p  C 

v 

irkjt 


kt  X2 
4o  C  t 

•  W  V 


(5.2) 


The  particular  solution  for  this  problem  can  then  be  obtained  by 
writing  T  =  Th  (x,t)F(t),  and  noting  that  conservation  of  energy 
requires 


/oo 

v.  T  dx  =  V 


(5.3) 


For  a  constant  heat  load,  the  wall  temperature  rise  obeys  the 
equation 


ATw  = 


2  %sft  _  r~ 

yfvp  C  kT  ^ 

v  w  V  T 


and  the  temperature  rises  at  a 
actual  rate  depends  on  the  material 
helps  select  among  various  candidate 


(5.4) 

rate  proportional  to  tH .  The 
properties,  in  a  way  which 
wall  materials. 


A  metal  which  is  a  good  candidate  for  the  bore  surface  of 
gun  barrel  should  have  a  large  heat  capacity,  a  large  thermal 
conductivity  and  a  high  melting  point.  The  merits  of  these 
characteristics  can  be  quantified  in  the  performance  factor,  f, 
defined  in  Equation  (5.4).  The  quantity  f  becomes  larger  for  a 
wall  material  which  is  more  resistant  to  heat  damage.  Table  1 
presents  experimental  values  of  key  transport  parameters  and 
computed  performance  data  for  various  barrel  conductors. 
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Material 

CU-AI2  O3 

Mo-Ti-Zr 

W-Ni-Cu 

W-Ni-Fe 

W-Ni-Fe-Cu 

Trade 

Name 

Glidcop 

AL-15 

TZM 

HD-17 

X-  5 

X-1.1 

Density  xlO3 

Kg/m3 

8 . 8 

9.0 

17 

18.4 

18 . 5 

Specific 

Heat  xlO3 

J/Kg-°  K 

0.50 

0.30 

0.16 

0.14 

0.15 

Yield  Strength 
MPa 

320 

380 

90 

120 

120 

Young ' s 

Modulus 

GPa 

120 

320 

280 

350 

350 

Thermal 

Cond.  xlOO 

W/m-° K 

3.2 

1.0 

0.94 

0.88 

0.88 

Electrical 
Conduct . 

%  IACS 

92 

32 

14 

13 

13 

Working  Temp. 

0  K 

1260 

1920 

1300 

1300 

1300 

Melting  Temp. 

1356 

1998 

1356 

1623 

1356 

Evap.  Temp. 

0  K 

2868 

5833 

6203 

6203 

6203 

Perf.  Factor  (f) 
xlO4  m2  -°  K/W  s 

0.30 

0.69 

0.70 

0.75 

0.72 

Work  Limit 
(Tw-To  )  /f 

3.2 

2 . 4 

1.4 

1.3 

1 . 4 

Melt  Limit 
(Tm-To  )  /f 

3.5 

2 . 5 

1.7 

1.8 

1 . 5 

Evap  Limit 
(Te -To  )  /f 

8 . 6 

8.0 

8 . 4 

HI 

8.2 

Table  1.  Experimental  and  computed  data  for  high-temperature  metal 
alloys  suitable  for  gun  barrels  and  rails. 
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The  first  two  materials  are  low-density  highly-conductive 
materials  appropriate  for  electromagnetic  gun  rails.  The  last 
three,  which  are  sintered  materials,  have  higher  density  and 
lower  electrical  conductivity,  and  may  be  appropriate  either  for 
a  segmented  rail  gun  or  for  a  liner  in  an  electrothermal  gun. 


The  material  properties  and  transport  coefficients  are 
needed  to  determine  the  behavior  of  the  gun  bore  under  the  heat 
load  and  stresses  delivered  to  it.  The  most  important  properties 
for  the  thermal  loads  are  the  working,  melting  and  vaporization 
temperatures . 


The  computed  performance  factor  (f)  determines  the  ability 
of  the  material  to  withstand  the  calculated  heat  load  from  the 
interior  plasma.  This  factor,  in  conjunction  with  the  three 
critical  temperatures,  gives  figures  of  merit  for  the  candidate 
materials  under  progressively  more  severe  conditions. 

The  figures  of  merit  are  summarized  in  the  Working,  Melting 
and  Vaporization  limit  factors  for  each  material.  They  show  that 
these  materials  may  soften  and  flow,  but  a  very  large  heat  input 
is  required  to  evaporate  them  and  contaminate  the  propelling 
plasma.  This  confirms  that  these  are  good  choices  for  gun  bore 
materials . 


Table  2  presents  comparable  data  for  selected  electrical 
insulators.  These  materials  are  necessary  in  a  rail  gun  in  order 
to  electrically  isolate  the  rails.  Unfortunately,  these 
electrical  insulators  are  also  exposed  to  the  propellant,  because 
they  and  the  rails  must  share  the  inner  surface  of  the  gun  bore. 
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Material 

B-N 

Si02 

AI2  O3 

Density  xlO3 

Kg/m3 

3 . 1 

3.2 

2 . 0 

2.2 

3 . 6 

Specific 

Heat  xlO3 

J/Kg-°  K 

1 . 3 

1 . 2 

2.3 

B 

1 . 2 

Flexural 

Strength 

MPa 

310 

700 

47 

80 

138 

Young ' s 

Modulus 

GPa 

410 

290 

80 

70 

280 

Compressive 

Strength 

MPa 

5160 

3500 

105 

1120 

345 

Thermal 

Cond .  xlOO 

W/m-°  K 

0.30 

0.18 

0.29 

0.05 

0.06 

Working  Temp. 

0  K 

2600 

1573 

2073 

1926 

2273 

Melting  Temp. 

0  K 

2773 

2173 

3255 

2001 

2326 

Evap.  Temp. 

0  K 

3073 

2173 

3255 

3223 

3253 

Perf.  Factor  (f) 
xlO4  m2  -°  K/W  s 

1.0 

1.4 

1.0 

3.2 

2.2 

Work  Limit 
(Tw-To  )  /f 

2 . 3 

0.9 

1.8 

0.5 

0.9 

Melt  Limit 
(Tm-To  )  /f 

2 . 5 

1.3 

3.0 

0.5 

0.9 

Evap  Limit 
(Te -To ) /f 

2 . 8 

1 .  3 

3.0 

0.9 

1 . 4 

Table  2.  Experimental  and  computed  data  for  high-temperature 
materials  for  rail-gun  electrical  insulators. 
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Comparison  of  the  insulator  temperature  limit  factors  with 
those  of  the  conductors  makes  it  clear  that  the  working  and 
melting  limits  are  in  the  same  range,  but  the  insulators  are  much 
more  easily  evaporated.  Other,  more  conventional  insulator 
materials  were  subjected  to  this  analysis  and  calculation,  but 
provided  much  worse  performance,  so  they  are  not  tabulated  here. 
The  key  results  for  them  are  summarized  in  figures  below. 

The  computed  performance  of  wall  material  under  a  heat  load 
is  illustrated  in  Figure  9.  The  "Scaled  Heat  Load"  is  (T-To)/f, 
and  can  be  compared  directly  with  the  limits  in  Table  1. 

This  omits  the  working-temperature  information,  but  compares  the 
melting  point  and  evaporation  point  data  with  the  measured 
experimental  mass  loss  of  gun  tube  liner  material.  The  onset  of 
mass  loss,  and  dramatic  increase  of  mass  loss  are  clearly 
correlated  to  the  calculated  melting  and  evaporation  points. 

The  analysis  and  simulation  of  this  process  has  suggested 
several  ways  to  reduce  the  barrel  erosion  problem.  It  is 
possible  to  add  materials  to  reduce  radiation  transport  in  the 
boundary  layer.  One  might  also  add  ablatable  materials  to  the 
wall  so  that  the  ablation,  when  it  occurs,  only  removes  the 
coating,  not  the  wall  material.  Finally,  since  so  much  of  the 
heat  transport  is  driven  by  the  fluid  temperature,  the 
performance  can  be  improved  most  directly  by  lowering  the  fluid 
temperature,  but  not  so  much  that  the  gun  performance  is 
degraded . 

The  20  mm  Electrothermal  gun  is  now  routinely  run  with  the 
hot  fluid  emanating  from  the  cartridge  cooled  by  mixing  with  a 
low- temperature  medium  which  is  placed  outside  the  capillary. 
The  mixed  lower-temperature  fluid  then  travels  down  the  barrel, 
separated  from  the  barrel  wall  by  a  very  thin  sheath  layer. 
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ABLATED  MASS  (mg) 


SCALED  HEAT  LOAD 


Figure  9.  Onset  of  Erosion  at  barrel  wall. 
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Figure  10  schematically  illustrates  the  experimental  method 
which  has  solved  the  bore  erosion  problem  for  a  wide  class  of 
Electrothermal  gun  applications.  The  solution  is  the  addition  of 
a  "mixing  chamber"  between  the  capillary,  which  produces  plasmas 
with  temperatures  of  30,000  0 K  to  40,000  °K,  and  the  gun  bore, 
which  requires  much  lower  temperatures. 

In  the  upper  diagram,  the  projectile  is  shown  separated  from 
the  capillary  by  a  chamber  filled  with  a  secondary  propellant, 
such  as  water  or  methanol.  When  electricity  is  fed  to  the 
capillary,  the  resulting  plasma  stream  must  pass  through  the 
mixing  chamber  on  its  way  to  the  bore. 

As  depicted  in  the  lower  figure,  the  fluid  mixes  with  the 
plasma,  flowing  with  it  into  the  bore  and  accelerating  the 
projectile.  In  practice,  the  mixing  fluid  contains  85-95%  of  the 
total  propellant  mass.  This  lowers  the  plasma  temperature  by 
more  than  an  order  of  magnitude,  reducing  it  to  temperatures 
comparable  with  those  in  conventional  guns. 

This  method  of  controlling  the  propellant  temperature  allows 
the  heat  transport  to  the  barrel  to  be  made  comparable  with  that 
in  conventional  guns,  while  the  choice  of  propellant  medium  (now 
primarily  the  mixing  medium) ,  still  produces  a  substantial 
improvement  in  gun  performance. 


BREECH  BLOCK 


PROJECTILE 


CONTAINED.  PROPELLING  FLUID 


FLUID 


Figure  10.  Solution  of  bore  erosion  in  Electrothermal  gun. 
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Electromagnetic  Results 


The  Electromagnetic  gun  may  be  thought  of  as  an 
Electrothermal  gun  in  which  the  capillary  runs  after  the 
projectile  and  pushes  directly  on  it  with  the  aid  of 
electromagnetic  forces. 

The  direct  action  of  the  electromagnetic  force  on  the 
projectile  is  a  clear  advantage.  It  is  not  necessary  to  send 
propellant  a. long  way  down  a  barrel  and  accept  propagation  losses 
before  the  propellant  can  deliver  energy  to  the  projectile. 

On  the  other  hand,  electromagnetic  radiation  ablates 
everything  present  in  the  capillary  of  an  ET  gun,  and  the 

armature  arc  in  the  EM  gun  is  ready  and  able  to  do  the  same 

ablating  in  the  bore  of  the  EM  gun. 

The  primary  heat  transport  in  this  case  is  radiation,  since 
the  hot  core  plasma  is  well  into  the  black-body  regime.  The  rate 
at  which  the  energy  reaches  the  wall  is  then  determined  by  the 
Rosseland  opacity.  The  best  possible  case  is  obtained  with 
Hydrogen,  for  which  the  Rosseland  opacity  is  shown  in  Figure  11. 
The  important  thing  to  note  is  that  this  quantity  has  a  very 

large  variation  in  its  values.  This  means  it  can  be  the  dominant 

process  in  one  region,  and  insignificant  nearby. 

For  the  more  restricted  range  of  temperatures  near  the  wall, 
we  have  developed  a  bilinear  fit  in  pressure  and  temperature,  so 
that  the  complete  Sesame  tables9  are  not  needed.  For  T  >  0.34  eV , 
the  fit  is 

loglO(LRoss)  =  0.35  -  P  *  0.35 

-(2.23  -  P  *  0.16  )  *  T 

+  (0.58  -  P  *  0.061  )  *  T*  *  2 

-(0.034-  P  *  0.0062)  *  T*  *  3 
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The  logarithm  is  continuous  at  T  =  0.34  eV,  and  is  proportional 
to  1/T  when  T  <  0.34  eV.  This  preserves  the  trends  of  the  Sesame 
tables  without  requiring  use  of  the  full  tables. 

Now  that  the  heat  load  is  known,  it  is  appropriate  to  return 
to  material  performance  data.  Figure  12  shows  the  mass  ablation 
data  for  candidate  insulator  materials.  In  addition  to  the 
materials  presented  in  Table  1,  popular  rail  insulator  materials 
such  as  Lexan,  Vespel,  Melamine  Glass  (G-9),  Epoxy  glass  (G-10) 
and  Polyimide  Glass  (PI-GL)  were  evaluated. 

Note  that  the  materials  of  Table  1  are  able  to  tolerate  heat 
load  for  some  time  before  melting  or  evaporating,  while  the 
composite  glasses  began  to  ablate  almost  immediately.  This  is 
why  it  is  important  to  examine  the  physical  properties  and 
evaluate  the  performance  factor.  It  is  a  useful  indicator  of  the 
suitability  of  materials  for  service  in  electrothermal  or 
electromagnetic  gun  bores. 

Figure  13  combines  the  results  for  rail  gun  conductors  and 
traditional  insulators  on  the  same  graph.  This  shows  that  the 
problem  of  insulator  ablation  is  much  more  severe  than  rail 
ablation,  when  conventional  insulating  materials  are  used. 

The  broken  line  at  18  milligrams  of  ablated  material  is  an 
estimated  lxmit  of  the  maximum  tolerable  amount  of  ablation  in  a 
small  rail  gun.  Clearly,  the  conventional  insulators  have 
exceeded  this  limit  before  any  high-performance  rails  have  shown 
any  significant  erosion. 
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The  only  insulators  that  have  any  significant  chance  of 
surviving  in  EM  service  are  Boron  Nitride,  Silicon  Nitride  and 
alpha  Silicon  Carbide.  These  materials  are  difficult  to  form 
into  appropriate  shapes  for  railgun  insulators,  but  there  are 
possible  designs  if  the  overall  gun  geometry  is  properly  chosen. 

The  fundamental  conclusion  is  that  the  conventional 
Electromagnetic  gun  will  operate  with  substantial  ablation  of 
both  rails  and  insulators  under  normal  conditions.  These 
materials  must  be  considered  as  the  constituents  of  the  armature 
after  the  projectile  has  begun  to  move,  and  also  as  parasitic 
mass  that  the  projectile  must  carry  with  it  as  it  accelerates 
down  the  bore . 
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Conclusions 


Radiation,  thermal  conduction  and  viscosity  and  turbulent 
convection  all  provide  important  transport  mechanisms  in  the 
plasma  flow  in  electric  guns.  The  flow  behind  a  projectile 
travels  at  a  velocity  near  the  speed  of  sound,  and  is  separated 
from  the  wall  by  a  very  thin  boundary  layer. 

The  heat  transport  through  this  boundary  layer  is  capable  of 
ablating  material  from  the  wall,  causing  erosion  which  damages 
the  gun  bore  and  reducing  the  performance  of  the  gun.  The  amount 
of  damage  this  transport  will  do  depends  on  the  temperature  and 
flow  momentum. 

In  the  Electrothermal  gun,  study  of  this  transport  problem 
has  led  to  a  redesign  of  the  basic  gun  system.  The  addition  of  a 
mixing  chamber  (Figure  10)  has  greatly  reduced  the  erosion 
problem.  The  choice  of  a  low-molecular-weight  propelling  fluid 
has  greatly  reduced  the  temperature  of  the  fluid  to  which  the 
wall  is  exposed  without  degrading  gun  performance. 

This  temperature  control  has  allowed  the  Electrothermal  gun 
to  operate  in  temperature  regimes  comparable  with  conventional 
guns.  However  the  ET  gun  performance  is  significantly  improved 
(Figure  5),  when  compared  to  conventional  gun  performance.  This 
is  because  of  its  ability  to  used  shaped  power  pulses  and  its  use 
of  low-molecular-weight  propelling  fluids  which  have  higher  sound 
speed  and  thus  more  efficient  coupling  from  the  breech  to  the 
projectile . 

In  the  Electromagnetic  (rail)  gun,  the  plasma  appears 
directly  behind  the  projectile.  Therefore  no  "mixing  chamber" 
design  is  possible.  The  insulators  ablate  as  a  result  of  the 
heating,  and  rails  may  also,  though  to  a  lesser  extent  (Table  1). 
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This  ablated  material  adds  to  the  accelerated  mass,  reducing  the 
efficiency  of  projectile  acceleration . 

Both  Electrothermal  and  Electromagnetic  guns  appear  capable 
of  reaching  higher  velocities  than  conventional  powder  guns. 
Both  have  demonstrated  this  in  the  laboratory.  However,  the 
relevant  measure  is  the  weaponizability  of  the  resulting  device. 

The  conventional  EM  gun  is  subject  to  a  wall  erosion  problem 
that  appears  unavoidable.  As  a  result,  military  operation  would 
require  replacement  of  the  barrel  after  a  few,  or  perhaps  up  to 
30  shots. 

The  recent  redesign  of  the  Electrothermal  gun  has  greatly 
reduced  barrel  erosion.  Present  designs  of  upgraded  ET  guns  use 
pressures  and  temperatures  comparable  with  present  120  mm  and 
155  mm  powder  guns,  and  shorter  acceleration  times.  Thus  barrel 
lifetimes  equal  or  longer  than  conventional  guns  are  predicted. 
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1.  INTRODUCTION 


A  simple  model  for  the  interior  gas  flow  inside  the  barrel  of  a  gun 
is  constructed.  The  approach  is  to  consider  the  case  of  a  projectile 
accelerated  at  constant  acceleration  due  to  a  constant  base  pressure,  and 
then  to  determine  the  gas  flow  that  must  be  supplied  at  the  gun  tube 
entrance  to  sustain  this  constant  acceleration.  The  effects  of  skin 
friction  between  the  gas  flow  and  the  bore,  as  well  as  thermal  transport 
from  the  gas  to  the  bore,  are  included.  Of  these,  skin  friction  has  the 
most  important  effect  on  gun  efficiency  due  to  its  contribution  to 
pressure  drops  along  the  gun  tube.  The  model  is  particularly  suited  for 
electrothermal  guns^  where  the  power  is  injected  electrically  into  the 
chamber  and  can  be  well  controlled  as  a  function  of  time. 

2.  ACCELERATING  PRESSURE  RAMP  MODEL 

We  seek  an  approximate  solution  in  which  the  projectile  velocity  Vp 
undergoes  uniform  acceleration,  g,  and  the  fluid  velocity  deviates  only 
slightly  from  spatial  uniformity  between  the  projectile  base  and  the  gun 
tube  entrance.  The  fluid  velocity  is  written  in  the  laboratory  frame, 

V  =  Vp(t)  +  v(x,t) 

Vp  =  gt  (1) 

V  gt2,t)  =0. 

It  is  convenient  to  view  the  pressure  profile  in  the  gun  in  a  frame  of 
reference  in  which  the  projectile  is  at  rest  at  the  origin  (Figure  1) 


Pressure  ramps  behind  a  uniformly  accelerating  projectile  (a)  in  the 
absence  of  skin  friction  and  (b)  with  skin  friction.  Note,  x  is  measured 
back  from  the  projectile  base  in  a  frame  at  rest  on  the  projectile  and  a 
uniform  gravitational  force  per  unit  mass,  g,  acts  on  the  gas. 


with  its  constant  base  pressure  Pproj-  In  this  frame  the  barrel  slides 
past  the  projectile  along  the  x  axis  at  velocity  gt  and  the  gas  experi¬ 
ences  a  gravitational  field 


A  P 
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where  A  =  ua2  is  the  gun  bore  cross-sectional  area. 


(2) 


2.1  Limit  of  Zero  Coupling  between  Gun  Tube  Mall  and  Gas  Flow 

In  this  simple  limit  the  only  force  on  the  gas  is  the  uniform 
gravitational  field  (2) so  that  the  pressure  ramp  in  Figure  la  is 
P  =  Ppr0jexP(x/h0) »  where  the  scale  height  hQ  =  c|/yg  with  c$  the  sound 
speed.  The  velocity  of  the  gas  in  the  laboratory  frame  is  V(t)  =  gt  with 
v  =  0  for  this  example,  and  provided  the  gun  tube  entrance  is  fed  with 
gas  of  pressure  P(x  =  gt2/2)  and  velocity  Vp(t)  the  self-accelerating 
ramp  of  gas  maintains  the  constant  pressure  Pproj  on  the  projectile  base. 


2.2  Model  Including  Skin  Friction  and  Wall  Losses 

(2) 

Gas-wall  interactions'  '  give  rise  to  important  effects  which  will 
now  be  included  approximately.  The  largest  effect  (assuming  no  ablation) 
is  the  skin  friction  force  per  unit  mass  acting  on  the  gas,  xV2/2d,  where 
the  friction  coefficient  X  is, 


X  = 


2  >°9U)  (if-)  +  ’-74 

'  w ' 


-2 


(3) 


and  6aw  is  the  roughness  scale  height  of  the  gun  bore  of  diameter  d  =  2a. 
This  friction  term  gives  rise  to  the  velocity  v(x,t)  in  equation  (1),  and 
both  x  and  v  are  treated  as  small  and  of  the  same  order,  i.e., 

X ~ v ~0( e) .  To  first  order  the  skin  friction  per  unit  mass  is  thus 
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XV*/2d.  It  behaves  like  a  slow  time  dependent  increase  in  the 
acceleration  factor  g  occurring  in  the  scale  height  formula.  Provided 
the  "exponential  atmosphere"  can  adjust  its  scale  height  h  rapidly  to 
these  slow  changes  in  g,  the  pressure  ramp  in  the  presence  of  skin 
friction  can  be  viewed  as  evolving  through  a  series  of  quasi-equilibrium 
states  given  by 


I  |£  =  _£  „p  =  g  + 

p  3X  y  3X 


XV2 


2d 


(4) 


so  that 


p(x,t)  pprojexP^  » 

1  cf  /  XV2  \ 

(h(t)  =  *  9*  =  9  v  +  29 §7 


(5) 


The  quasi-equilibrium  model  (4)  applies  provided  the  scale  height  h 
becomes  shorter  due  to  skin  friction  with  a  speed  satisfying  dh/dt  <  c$ 
so  that  the  pressure  ramp  (5)  can  adjust  via  local  sound  wave  propagation. 
This  condition  becomes  t  <  yd/\cs  after  which  time  Vp  increases  to  the 
point  where  the  ramp  can  no  longer  completely  adjust. 

In  order  to  evolve  the  series  of  ramps  (5)  the  gas  must  have  a 
velocity  v  (in  the  frame  of  Figure  lb)  given  by 


v(x’t}  =  (ll)n~  “("/)  x  * 


(6) 


The  entrance  pressure  P£nter  then  evolves  along  the  dotted  path  shown  in 
Figure  lb. 
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Heat  loss  from  the  gas  to  the  bore  is  also  occurring  but  has  a  less 
important  effect  on  the  scale  height  h.  It  will  be  accounted  for  later 
in  the  energy  balance. 

3.  ENERGY  INPUT  AND  EFFICIENCY  OF  THE  GUN 

We  now  view  the  gun  in  the  laboratory  frame  (Figure  2).  Gas  flowing 
into  the  gun  tube  entrance  must  carry  sufficient  energy  to  "fill  out"  the 
series  of  self-accelerating  quasi-equilibrium  ramps  (5).  We  assume  this 
input  power  continues  for  a  pulse  time  t  after  which  it  tapers  off.  The 
total  energy  that  resides  in  the  gun  at  a  time  t  <  t  can  thus  be  written 
using  (5)  and  (6)  and  integrating  from  x  =  0  to  xp  =  V2/2g  for  the  barrel 
contributions. 


W. 


A  Pr 

input  (V 

(gas  internal  energy) 


^  (*v> 


(*"•'"  -<)(<■  ¥) 


2xx  x/h 

e  P 


(7) 


(gas  kinetic  energy) 


+  a  p  y  A  xeffPEnter^  +  nm 

f  A  HProjxp  +  - (7TT) - 

(projectile  KE)  (energy  in  chamber) (loss  to  wall) 
where  the  effective  chamber  length  xg^  is  defined  by 

A  xg^  =  ira2xe^  =  chamber  volume,  (8) 

and  Pproj  =  YPproj/c|.  The  entrance  pressure  into  the  gun  tube  is  from  (5), 
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Fiqure  2 


Parameters  for  an  accelerating  gas  ramp  in  a  gun  tube. 


yVP 

PEnter  PProj  exp  ^ 


(>*S) 


and  in  the  laboratory  frame  (Figure  2)  the  gas  flow  follows  from  (6)  as. 


W  ■  vp +  25&  •  001 

and  pressure  drops  between  the  chamber  interior  and  gun  tube  entrance  are 
neglected  in  the  chamber  energy.  It  is  also  interesting  to  note  that  if 
we  use  VEnter  and  PEnter  to  calculate  a  power  flow  into  the  gun  tube  and 
integrate  it  from  0  to  t,  a  result  slightly  different  from  the  barrel  part 
of  (7)  is  obtained.  This  is  because  the  quasi-equilibrium  ramps  are  not 
exact  solutions  to  the  flow  equations. 


To  evaluate  thermal  loss  to  the  tube  wall  we  use  the  convective 
(21 

:nnv*f  fnrmnl  a  '  ' 


transport  formula. 


q  =  *PV 


ergs/cm2/sec  , 


which  gives 


t  x„(f) 


Q(t)  •  /dt/ dx  q  T,  KrV  Pproj  je  -  1  - 


0  0 
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where  we  assumed  (T  -  TWa-^)/T  «  1  (T  =  gas  temperature).  The  function  $ 
appearing  in  (11)  and  (12)  is  discussed  in  Ref.  2,  Eq.  (3.19),  and  is  about 
1.5  for  neutral  gas  and  smooth  bore  gun  tubes.  Heat  transfer  to  the 
chamber  wall  will  be  neglected. 


3.1  Choked  Flow  Limit  to  V 


Enter 


The  maximum  velocity  with  which  gas  can  flow  into  the  barrel  from  a 
chamber  of  chambrage  >  1  is  the  sound  speed.  The  steepening  pressure 
ramps  can  thus  continue  to  be  fed  at  the  entrance  provided  from  (10), 


mc<) 


i.e.,  the  pulse  time  t  must  satisfy 


T  <  T  = 

"  choke  g  \cs) 


choke 


(13) 


where 


( 


$ 


s  /  choke 


1  -»r 


1/3 


(14) 


with 


=  ds_ 


(15) 


AC, 


The  quantity  (V  /cs)choke  is  graphed  as  a  function  of  y  in  Figure  3. 
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Projectile  velocity  V  divided  by  the  sound  speed  c  at  the  time  in  its 

P  ^ 

acceleration  that  gas  flows  into  the  gun  tube  entrance  at  the  sound  speed. 
Note  that  a  small  acceleration,  g,  and  a  large  skin  friction  coefficient 
causes  the  projectile  velocity  to  reach  only  a  small  fraction  of  c$  when 
the  choked  flow  limit  c  is  encountered  at  the  gun  tube  entrance. 


3.2  Rarefaction  Time,  Projectile  Energy,  and  Efficiency 

Projectile  acceleration  continues  beyond  the  power  input  time  t. 

Assuming  the  input  power  is  reduced  rapidly  after  t,  a  rarefaction 
travels  from  the  chamber  down  the  gun  tube  to  the  projectile  base  where  a 
deeper  rarefaction  is  reflected.  In  order  to  launch  a  strong  rarefaction 
into  the  barrel  (so  that  no  further  appreciable  acceleration  occurs  after 
its  arrival  at  the  projectile  base)  we  first  note  that  it  takes  a  time 
xe^/2V£nter  for  the  chamber  to  drop  its  pressure  by  0.5  after  pulse 
shutdown.  Since  the  projectile  has  advanced  to  a  distance 

O 

g(x  +  xeff/2VEnter)  /2  by  this  time,  the  total  time  xRare  for  the  rarefaction 
to  reach  the  projectile  is 


Rare 


(t  + 


Vff\ 


eff 


2V 


Enter 


2  V 


Enter 


/  V  \ 

2\cs  ) 


xeff(y  4  cs> 
2cs  ^Enter.t 


06) 


where  the  added  subscript  x  means  quantities  are  evaluated  at  time  x. 


The  total  energy  acquired  by  the  projectile  is  thus 


W  =  —  M  V2 
wProj  2  nProj  VF  * 


(17) 


where  the  projectile's  final  velocity  is 


VF  =  9(x  +  TRare}  ’ 


(18) 


which  using  (10)  and  defining 


6  =  ^ 


(19) 


gives 


9 


2c 


1  +  6 


6  + 


Xcs 

2gd 


B' 


(20) 


and  we  note  6  has  a  maximum  given  by  (14)  for  pulses  of  maximum  length 
Tchoke* 

The  gun  efficiency  can  now  be  written  as 


w 

=  wProj _ 

Winput^  +  W1ate 


(21) 


where  WinpUt;  is  given  by  (7)  and  W-jate  is  the  energy  input  to  the  gun 
occurring  after  the  peak  power  input  time  x. 

4.  SUMMARY  OF  FORMULAS 

It  is  useful  to  define  some  dimensionless  quantities 


e  = 


eff 

V 


(jj  = 


AxC< 

“2T 


=  _2l  =  51 
cs  cs 


(22) 


The  efficiency  (21)  then  becomes. 


_  0.5(y-1)f  VF 


y(1  ±  ^b) 
(l  ♦ 


(23) 


W 

inputs 


A  = 


2.10_3<t>  fxCsi 


1  +  7  ^(T-l  )H  -  £  ♦  2B*]  +  obyO+*b)  +  |  X 


exp  ( 1  +ipB )J  +  y ( Y-l  )Bz(1+i|/B)  -  1  -  y  y62(y-1)(1  - 


(24) 


2.1 0-3<t>  TCs 


[tttVbT  +  ?  *'*]} 
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and  the  projectile  final  velocity  Vp  is  given  by 


+ 


e_ 

2 


1  +  B 

1  +  fty 


and  the  formulas  only  apply  for  pulses  of  length 


(25) 


2xeff 

Bcg  T  "  Tchoke 


- 

^  \  cs  /  choke 


(26) 


given  by  (14),  (15),  and  Figure  3.  The  lower  limit  on  t  in  (26)  exists 
because  the  barrel  hydrodynamics  have  been  treated  more  accurately  than 
the  chamber  hydrodynamics  in  this  model.  The  condition  is  that  the  barrel 
volume  V^tA/2  should  exceed  the  chamber  volume  xeppA  for  the  model  to 
apply,  e.g.,  by  a  factor  of  3  or  more. 

Other  quantities  of  interest  are  the  pressure  ratio 


PEnter^ 

p 

Proj 


(27) 


and  the  total  mass  of  working  fluid  in  the  gun 


'FI 


uid  * 


PT 


I  P.(x). 


dx  +  Ax 


yP 


Enter, t 


eff  \  c2 


(28) 


Proj 


[t6B  +  TTT*bt]  exp  (?  (1+Wi))  ■ 


Mfi  ..  fixes  the  minimum  chamber  size  needed  to  contain  this  mass  of 
r  luid 
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fluid  at  its  original  liquid  or  solid  density. 


Finally,  consider  the  power  input  to  the  gun.  At  the  end  of  the 
pulse  time  t  we  assumed  the  fluid  in  both  the  barrel  and  the  chamber  to 
have  the  same  y.  However,  at  early  times  the  fluid  in  the  chamber  may 
have  a  different  value,  say  y  ,  due  to  EOS  effects  such  as  molecular 
covolume.  If  we  assume  that  the  very  early  part  of  the  power  input  pulse 
raises  the  chamber  volume  to  pressure  Pproj  before  the  projectile  advances 
appreciably  (say  not  more  than  about  one  bore  diameter),  then  the  above 
formulas  are  not  sensitive  to  the  fill-up  pulse  shape.  The  energy  required 
in  this  early  pulse  is  ^xeffPproj/(Y0"^ )  anc*  should  be  delivered  in  a  time 
(Figure  4). 


After  this  the  power  input  required  is  given  by  Winput  from  (7)  until  the 
end  of  the  effective  pulse  time  t.  Energy  flowing  in,  Wlate,  for  t  >  t 
•  has  little  effect  on  the  projectile  due  to  the  rarefaction  catch-up  and 

reflection.  This  late  energy  simply  reduces  the  efficiency  (23)  or  (21). 


4.1  Gun  Tube  Entrance  Heating 

The  gun  tube  entrance  has  a  wall  temperature  Tw  and  receives  heat 
from  the  gas  flow  of  temperature  T.  Using  Eq.  (4.1)  in  Reference  2,  this 
increases  to  Tw(t)  °K  given  by 


v°> 


dt 


(T-t) 


104  y 

1/2  vm 


Pkb^km/s/Enter 


t  -  v*) 
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Tfin 


Tchoke 


Figure  4 


Input  power  to  a  gun  showing  the  three  components,  namely  (i)  pressure 
fill-up  of  the  chamber,  (ii)  the  constant  acceleration  gas  ramp,  and 
(in)  the  late  wasted  energy  input  the  beginning  of  which  launches  a 
rarefaction  into  the  barrel.  Note  W-|  te  includes  all  the  energy  input 
for  times  t  >  smaller  of  x  or  x^^  and  the  case  x  <  x^^  is  shown. 


in  the  pulse  time  t.  We  evaluate  this  in  the  same  approximate  manner  as 
Reference  2  and  find*using  (9)  and  (10), 


Two  +  T(T/Theat>,/2 


where 


0.525  (y-1) 

(T(°K)/3,500) 

• 

Theat 

ccy<i>B  C$(km/s) 

PProj (kbar,)exp[If2;  (!+♦&) 

2 

seconds. 


(30) 


(31) 


This  heating  is  severe  and  leads  to  the  desirability  of  wall  coatings  at 
large  bore  sizes.  High  sound  speeds,  low  pressures,  low  gas  temperatures, 
high  y  and  small  a  are  also  desirable  (a  =  2/(irpwc  1^)^  in  units  of 
cm2  °K  sec^2/Joules) .  But  low  pressures  cause  low  accelerations  and  the 
limits  imposed  by  skin  friction  and  choked  flow  (eq.  14)  lower  the  gun 
performance  and  must  be  involved  in  the  trade-offs. 


5.  NUMERICAL  DISCUSSION 

The  sound  speed  cs  and  specific  heat  ratio  y  are  important  quantities 
for  the  propelling  gas.  The  sound  speed  strongly  determines  the  scale 
height  h  of  the  gas  ramp  (Eq.  5)  which  in  turn  limits  the  practical 
acceleration  length  for  a  given  accelerating  pressure  pproj »  the  limit 
being  set  by  the  ability  of  the  chamber  to  contain  the  entrance  pressure 
and  by  the  increased  bore  entrance  heating  (Eq.  30)  that  occurs  if 
PEnter  becomes  very  lar9e-  The  skin  friction  dependence  of  h  (Eq.  5) 
also  sets  velocity  limits  for  low  mass  projectiles  accelerated  over 
long  lengths  by  low  propelling  pressures. 

*In  Equations  4.5  and  4.6  of  Ref.  1  set  t=i  and  t=t.  to  align  these 
results,  but  note  the  result  (31)  accounts  for  the  ba^fil  pressure  ramp. 
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Figure  5  shows  a  plot  of  efficiency  as  a  function  of  y  'For  a 
specific  example.  The  model  we  have  constructed  assumes  that  both 
y  and  cs  are  given  constants.  However,  they  must  be  chosen  to  be 
consistent  with  the  appropriate  equation  of  state  (EOS)  for  the  working 
fluid  and  for  the  pressures  and  temperatures  involved  in  the  gun.  Since 
the  highest  fluid  pressure  and  energy  density  occur  at  the  gun  tube 
entrance  at  the  end  of  the  pulse,  and  this  also  controls  the  chamber 
energy  and  choking  limit,  the  values  of  y,  cs,  evaluated  at  pEnter(0 
provide  a  reasonable  approximation  for  this  choice.  An  iteration  pro¬ 
cedure  could  be  used  between  the  model  equations  (23)  -  (26)  and  EOS 
to  select  y,  c$  consistently  with  pEnter  and  T.  In  the  rest  of  this 
discussion  we  give  an  example,  but  use  a  simpler  first  step  approxima¬ 
tion  to  such  an  iteration  process. 

Figures  6  and  7  show  the  sound  speed  c$  and  (y-1)  for  water  (super¬ 
heated  steam)  graphed  as  functions  of  temperature  for  several  pressures^. 
The  variation  of  (y-1)  with  P  is  relatively  more  important  as  can  be 
seen  by  noting  that  for  a  fixed  temperature,  e.g.,  T  =  3,273°K  (0.28  eV) 

(y-1)  «  0.169  +  0.021  Pkb  ,  (32) 

cs(km/sec)  «  1.356  +  0.074  Pkb  . 

To  include  this  approximately  in  the  expressions  for  e,  W.  £,  etc.,  we 
will  replace  y  by  'Y£nter  in  the  sensitive  parts  of  the  expressions  (on 
the  left  side  of  the  equations  below)  and  set  y  =  y,  cg  =  c$  on  the  right 
where  y,  £'  are  estimated  means  values  for  the  insensitive  parts  of  e. 
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lYEnter"^ 


0.5  y ( 1  +1^3 ) 


w 

input 


^input^YEnter"^  1  -2  / 
MProj  "  2  '*  U 


YEnter”^ 


PEnter  *  PProj  exp  i  *2  ^  +  <*> 


g  =  ^A 


All  quantities  on  the  right  are  functions  of, 


a’  T’  xeff ’  *»  X’  Y’  V  Mn _ . 


in  which  y,  c  ,  are  estimated  values  and  on  the  left  of  (33)  and  (34), 


YEnter  Y^PEnter^ 


where  the  example  (32)  for  3,273°K  steam  will  be  used.  It  should  also 
be  noted  that  for  pulse  times 

T>,ch0^yv3|[(,t|,)-1]’/3-[(1^r-'],/3i. 


y  =  dg/xc; 
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the  energy  in  that  part  of  the  pulse  W(t  >  *ch0ke)  is  not  effectively 
used  for  projectile  acceleration  and  should  be  included  in  W-|ate.  For 
these  cases  we  simply  set 

9Tchoke 


3  =  3 


choke 


(40) 


throughout  the  formulas. 

Using  the  above  procedure,  we  have  graphed  several  quantities  of 
interest.  Figure  8  shows  the  projectile  velocity  as  a  function  of  projec¬ 
tile  mass  for  a  500  ysec  input  pulse  time  for  the  20  mm  gun  parameters 
shown.  All  of  the  figures  9  to  15  give  results  for  the  optimum  pulse 
length  t  =  f°r  a  9un  with  parameters  listed  in  Figure  9.  These 

examples  are  all  for  superheated  steam  at  3,273°K. 

5.1  Model  Limitations  and  Chamber  Hydrodynamics 

The  model  gives  attention  to  the  pressure  ramp  in  the  barrel  and  the 
entrance  flow,  with  only  rough  accounting  for  the  chamber  hydrodynamics. 
Thus  the  model  has  useful  accuracy  for  guns  with  small  enough  chambers  so 
that  the  last  term  in  Eq.  20  is  small.  However,  a  more  exact  treatment  of 
the  chamber  hydrodynamics  would  require  one  to  address  a  specific  chamber 
since  gun  chambers  are  varied  in  geometry  and  internal  structure.  This 
is  particularly  the  case  for  electrothermal  guns  that  are  designed  to 
inject  these  constant  acceleration  ramps  into  the  gun  tube.  Note  also 
that  an  ideal  power  input  pulse  would  be  sharply  cut  off  at  t  as  in 
Figure  15.  However,  realistic  input  power  pulses  would  not  cut  off  so 
sharply  and  continuing  power  input  would  slightly  delay  the  rarefaction. 
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The  50%  chamber  pressure  drop  time  used  in  Eq.  16  is  an  estimate  for  all 
these  effects. 

A  second  area  where  the  model  could  be  improved  would  be  to  allow  the 
parameters  y  and  cg  in  the  scale  height  formula  (5)  to  vary  slowly  with  x. 
This  would  allow  one  to  account  for  the  pressure  dependence  of  these 
quantities  along  the  pressure  ramp  via  the  equation  of  state,  and  to  also 
allow  for  a  variation  in  the  gas  temperature  T(t)  injected  by  more  general 
pulse  structures. 


Equation  of  State  data  for  water 
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Projectile  velocity  as  a  function  of  projectile  mass  for  several  base  pressures  for  input 
power  pulses  of  length  t  =  t  .  .  .  The  remaining  figures  10  -  15  give  other  quantities 
for  this  same  example. 


Peak 

entrance 

pressure 

kbars 


entrance  pressure  as  a  function  of  projectile  mass  for  several  projectile  base  pressures 
ulses  of  length  x  =  xrV,_.,0  with  gun  parameters  the  same  as  in  Figure  9. 


Input 

Energy 

kilojoules 
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Input  Energy,  WjnpU^»  as  a  function  of  projectile  mass  for  several  projecti 
pressures  and  pulses  of  length  x  =  x  .  .  with  gun  parameters  the  same  as  1 
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